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ABSTRACT 

We present results from a mosaic of nine Chandra observations of M86 and the surrounding field. 
We detect three main diffuse components: the Virgo ICM at ~2.4 keV, the extended halo of M86 at 
~1.2 keV, and the cooler central and stripped gas of M86 at ~0.8 keV. The most striking feature is a 
long tail of emission, which consists of a plume ^ 4'north of M86 and two main extensions emanating 
from the plume. Based on the morphology and temperature structure of the tail, we conclude that 
it is formed by ram pressure stripping of M86 as it falls into the Virgo cluster and interacts with the 
Virgo ICM, in agreement with earlier work. The tail is 150 kpc in projection, and a simple estimate 
gives a lower limit on the true length of the tail of 380 kpc, making this the longest ram pressure 
stripped tail presently known. The total gas mass in the plume (~ 7 x 10® Mq) and tail (~ 1 x 10^ M©) 
is about three times that in the core of M86, which supports the scenario where most of the gas was 
stripped rapidly and recently. The projected position of the plume can be understood if M86 has an 
aspherical potential, as suggested by optical isophotes. Ram pressure stripping from an aspherical 
potential can also explain the split "double tails" seen in M86 and in other Virgo cluster galaxies 
in the field. The large line-of-sight velocity of M86 (1550 km s~"'^with respect to M87), its position 
relative to the Virgo cluster, and the orientation of the tail tightly constrain its orbital parameters. 
The data are inconsistent with a radial orbit, and imply inner and outer turning radii of w 300 kpc 
and To > 8.8 Mpc, indicating that M86 is, at best, only weakly bound to the Virgo cluster. 
Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Virgo) — magnetic 
fields — X-rays: galaxies — galaxies: individual (NGC4406) — galaxies: individual 
(M86) 
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1. INTRODUCTION 

M86 (NGC 4406) is a bright eUiptical (E3/S0) galaxy 
in the Virgo cluster of galaxies. It is the dominant 
member of one of the larger subgroups within Virgo 
(Binggeli et al. 1993; Bohringer et al. 1994; Schindler et 
al. 1999). Its line-of-sight velocity relative to M87, the 
dominant member of the Virgo cluster, is -1550 km s~^, 
much higher than the average cluster velocity dispersion 
(Smith et al. 2000). The X-ray surface brightness dis- 
tribution is unusual, with a large "plume" extending to 
the northwest from M86, which was first noticed as part 
of a survey of Virgo cluster galaxies undertaken with the 
Einstein Observatory (Forman et al. 1979). The galaxy 
has an optical asymmetry that extends in a direction 
similar to the direction of the plume (Nulsen & Carter 
1987; Mihos et al. 2005). Several authors have inter- 
preted this plume as arising from ram pressure stripping 
due to strong interactions with the Virgo ICM (Forman 
et al. 1979; Fabian et al. 1980; Takeda et al. 1984; Knapp 
et al. 1989; Bregman & Roberts 1990; White et al. 1991; 
Rangarajan et al. 1995). Elmegreen et al. (2000) find 
dust streamers in the core of M86 that connect to the 
nucleated dwarf galaxy VCC 882, and suggest a recent 
interaction between the pair that may have contributed 
to the asymmetry of M86's optical isophotes. 

At X-ray energies, M86 has been observed by Einstein 
(Forman et al. 1979; White et al. 1991), Ginga (Takano 
et al. 1989), EXOSAT (Edge 1990), ASCA (Matsushita 
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et al. 1994), and ROSAT (Bohringer et al. 1994; Ran- 
garajan et al. 1995). Recently, Finoguenov et al. (2004) 
presented XMM-Newton observations of M86 and con- 
cluded that its unusual morphology is due to an interac- 
tion with an X-ray filament rather than the Virgo ICM. 
However, their conclusions were based on the large sep- 
aration between M86 and M87 along the line of sight of 
2.4 ± 1.4 Mpc reported by Neilsen & Tsvetanov (2000), 
which is inconsistent with the more recent result from 
Mei et al. (2007) who find a separation of 0.4 ±0.8 Mpc. 

We report here on a mosaic of nine Chandra observa- 
tions of M86 and the surrounding field, totaling 240 ksec 
of exposure. The observations and data reduction tech- 
niques are described in §[2l The X-ray image is presented 
in § [21 and results on temperature and abundance struc- 
ture from spectral analysis are given in § [4l In § [5l we 
discuss some of the more interesting features of M86, and 
place constraints on its orbit. In particular, we give new 
results on the extent of the stripped tail, and measure 
the density profile of an X-ray brightness edge seen to 
the southeast. Our results are summarized in §[6l 

We assume a distance to the Virgo cluster of 16 Mpc 
throughout, consistent with the latest results from Mei et 
al. (2007), which gives a scale of 0.08 kpc/" for = 0.3, 
fl\ = 0.7, and Hq = 70 km s^^ Mpc^^. All error ranges 
are 90% confidence intervals, unless otherwise stated. 

2. OBSERVATIONS AND DATA REDUCTION 

Table [T] summarizes the nine Chandra observations of 
M86, which completely cover the center of the galaxy, as 
well as the plume and long tail. All data were reprocessed 
from the level 1 events files using the latest calibration 
files (as of CIA03.3). CTI and time-dependent gain 



2 



RANDALL ET AL. 



corrections were applied where applicable. LC_CLEAN 
was used to remove background flares'^. After periods 
with obvious flares were removed, the mean rate was cal- 
culated for some observed quiescent time interval, and 
the program was re-run, forcing the global mean to equal 
this quiescent rate. For observations without obvious 
flares, time bins that were not within 3cr of the mean 
were discarded. The final cleaned exposure times are 
given in column (5) of Tabic [T] 

The emission from M86 and the surrounding Virgo 
cluster fills the image field of view for each observa- 
tion. We therefore used the standard CALDB'' blank 
sky background files appropriate for each observation, 
normalized to our observations in the 10-12 keV energy 
band. We combined background maps for each pointing 
separately, and generated exposure maps to account for 
the pointing offsets and the different CCD responses. To 
generate exposure maps, we assumed a MEKAL model 
with kT = 1 keV, Galactic absorption, and abundance 
of 30% solar at a redshift z = 0, which is consistent with 
typical results from detailed spectral fits (see § . 

3. THE X-RAY IMAGE 

The exposure corrected, background subtracted, 
smoothed mosaic image is shown in Figure [T] Several 
patches of diffuse emission are seen in M86, as is an im- 
pressive tail extending to the northwest. To enhance 
the visibility of the diffuse emission we created an im- 
age with bright point sources removed. For each Chan- 
dra pointing, regions containing point sources were "filled 
in" using a Poisson distribution whose mean was equal to 
that of a local annular background region. The resulting 
smoothed image (Figure [2]) shows: 

• a long tail extending to the NW, first detected by 
Forman et al. (1979). The tail is bifurcated at its 
base, with the brighter (northern) extension at a 
position angle of 305° (measured from north to 
east) and the shorter extension at 285°. At 18.4' 
(88.4 kpc) from the center of the bright emission 
at the base of the tail (or 21.1' (101.5 kpc) from 
the center of M86), the tail turns directly north 
for 6.8' (32.5 kpc). It then curves to the west at 
26.7' (128.1 kpc) from M86, and finally is no longer 
detected at 30.6' (150 kpc). There is evidence for 
a fainter, parallel tail along the northeastern edge 
of the tip of the bright tail, discussed below. 

• a large plume of emission, directly north of M86, 
at the root of the long tail (noted previously by 
Forman et al. 1979). The tail emanates from this 
clump, rather than from M86 itself as one would 
expect for a continuous stripping process. 

• a sharp boundary along the northeastern edge of 
the tail, and an apparent "void" in the diffuse emis- 
sion just north of the tail. This is in stark contrast 
to the southwestern edge of the tail, where the sur- 
face brightness falls off smoothly. Rangarajan et al. 
(1995) first noted this void, and suggested that it 
is an evacuated Mach cone left from M86's passage 
through the outer regions of the Virgo ICM. 
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• an extended halo of X-ray emission associated with 
M86 with an edge 13' (62.4 kpc) to the southeast. 
This edge also is visible in ROSAT observations 
(see Figure [3]). 

• an asymmetry in the diffuse emission near the core 
of M86, with an extension to the south. Rangara- 
jan et al. (1995) suggest that this is a high pressure 
region, formed as the leading edge of M86 plows 
into the Virgo ICM. 

• several nearby galaxies. The elliptical galaxy M84, 
west of M86, shows complicated structure in the 
diffuse emission in the core related to the ra- 
dio outburst, and a tail extending to the south 
(Finoguenov & Jones 2002). Also visible are 
NGC 4388 (south-southwest of M86; Beckmann et 
al. 2004) and NGC 4438 (east of M86; Machacek 
et al. 2004; VoUmer et al. 2005). Both also have 
diffuse tails, which likely indicate their projected 
direction of motion with respect to the Virgo ICM. 
Most tails split into two main streams. A small 
patch of faint emission also can be seen in the re- 
gion of the disk galaxy NGC 4402 (Growl et al. 
2005). 

For comparison, the (heavily binned) Chandra mosaic 
image is plotted alongside the ROSAT and DSS images 
in Figure [3] The ROSAT image clearly shows the distri- 
bution of diffuse emission on larger scales. In particular, 
the contribution from the M87 halo is seen in the south- 
east, as is a large halo of diffuse emission around M86 
itself, outside the brightest central regions. These com- 
ponents are individually analyzed in our detailed spectral 
fits (see ^ 14. 2p . Results from these fits indicate that the 
extended M86 halo is composed of group gas associated 
with M86, at a temperature of kT « 1.2 keV. 

A close up view of the stripped tail, with scaling and 
binning chosen to show its structure more clearly, is pre- 
sented in Figure U] The twisting and curving of the tail 
at the faint tip may be caused by turbulence or similar 
"weather" in the Virgo ICM. The tail is clearly bifur- 
cated at its base. Additionally, a fainter tail, split from 
the first, is barely visible along the northern edge of the 
faint end of the main tail. To test the significance of 
this feature, we extracted the total count rate in evenly 
spaced bins across the width of the box shown in Fig- 
ure H) Figure [H] shows two distinct peaks. If we take an 
average count rate from the highest three bins in each 
peak and compare it with the average count rate from 
the lowest two bins in the gap, we find that the gap is 
significant at about 2.5a. As described below in § 15.41 
these double streams, as seen clearly at the base of the 
tail and possibly at the faint end of the tail, as well as in 
other Virgo cluster galaxies in the field, may be explained 
as stripping from an inclined aspherical potential. 

4. SPECTRAL ANALYSIS 

In summary, the X-ray images show three diffuse emis- 
sion components: the Virgo ICM (centered on M87), the 
extended M86 halo, and the core and tail of M86. We 
generate a temperature map as a guide for detailed spec- 
tral fitting to disentangle the various components. We 
assume a galactic absorption of Nh — 2.62 x 10^" cm~^ 
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throughout. AUowing the absorption to vary did not sig- 
nificantly improve any of the spectral fits. 

4.1. Temperature Map 

The temperature map was derived using a method em- 
ployed by O 'Sullivan et al. (2005) and Maughan et al. 
(2006). For each temperature map pixel (19.7"/pix), we 
extracted a spectrum from a circular region containing 
1000 net counts (after subtracting the blank sky back- 
ground and a component to account for emission from 
M87, see § 14. 2|) . up to a maximum radius of 3.3' (on 
the order of the chip size). This large maximum radius 
was used to obtain sufficient source counts to extend the 
temperature map to faint regions of the tail. The result- 
ing spectrum was fit in the 0.6 - 2.0 keV range with an 
absorbed APEC model using XSPEC. Data from back- 
side- (BI) and front-side-illuminated (FI) CCDs were 
treated as separate data groups for each observation. For 
each fit, the abundance was fixed at 26% solar (consistent 
with most of the detailed spectral fits; see Table [2|). The 
resulting temperature map is shown in Figure[6] (middle) . 
Regions that did not have at least 1000 source counts 
within the maximum 3.3' radius were excluded from the 
map, as were regions in the southeast that were dom- 
inated by Virgo cluster emission. For comparison, we 
also show the tessellated temperature map in Figure [HI 
Each bin was fit using only counts from that area, so the 
extraction regions are well-defined. The bins were gen- 
erated using the algorithm provided by Diehl & Statler 
(2006), which is a generalization of Cappellari & Copin's 
(2003) Voronoi binning algorithm, and requiring roughly 
1100 net counts per bin. Each bin was fit with a single 
APEC model, without a component to model the Virgo 
ICM background, in contrast to the "smoothed" tem- 
perature map. The size of the bins roughly indicate the 
size of the extraction regions for the smoothed temper- 
ature map pixels in the same area. The tessellated and 
smoothed temperature maps are in good agreement, ex- 
cept in the regions far from M86 where the contribution 
from the Virgo cluster is more important. 

In these temperature maps, the remarkably sharp 
boundary between the cool galaxy gas and that of the 
cluster on the northern edge of the tail is apparent (Fig- 
ure O. At least four distinct cold clumps can be seen: 
one centered on M86, one 1.41' (6.8 kpc) south of M86, 
one elongated clump 2.45' (11.7 kpc) east of M86, and 
a larger cool plume 3.63' (17.4 kpc) to the north. Each 
clump corresponds to a bright peak in the X-ray im- 
age. Additionally, the cooler tail extends to the north- 
west. To the southeast, in the region of the bright- 
ness edge discussed in § 15.21 the gas is somewhat hotter 
(~ 0.25 keV) in the southern half of this feature than 
in the northeastern half at the same distance from M86. 
This temperature difference is detected in spectra from 
non-overlapping regions in these two areas, indicating 
that the difference is not an artifact of the technique 
used to generate the temperature map. More detailed 
fits are discussed in § 14.21 

4.2. Detailed Spectra 

Based on the derived temperature map and the X-ray 
image, we defined 23 regions for detailed spectral anal- 
ysis (see Figure [7]) . To include data from a particular 
observation in the spectral analysis, we required that 



it contribute at least 200 counts (each data group was 
binned such that there were at least 40 counts per bin). 
A summary of the regions and spectral model for each 
region are given in Tabled Rl - R4 are centered on the 
individual clumps of emission near M86, while R5 -RIO 
trace the long tail NW of M86 (hereafter regions will be 
labeled as 'R#', with the number defined by column 1 
of Table [2]). Rll roughly corresponds to the "shocked 
region" indicated by Finoguenov et al. (2004; see their 
figure 7). There are several regions to either side of the 
M86 tail. R17, R18, R19, & R23 measure the large scale 
diffuse components from M86 and M87, outside the cen- 
tral region where the cool gas is observed (note that R23 
is a combination of five disjoint regions, see Figure [7]). 
Finally, R14 - R16 contain M84 and its tail of diffuse 
emission. Each region was initially fit with an absorbed 
APEC model over the 0.6-5.0 keV range, with an ad- 
ditional APEC component to model emission from the 
Virgo ICM as in § I4.H with the abundance allowed to 
vary. Additional components were added, as needed, to 
improve the fits (see Table [21 For example, a single tem- 
perature model gave large residuals at high energies when 
fit to the regions corresponding to the cores of M86 and 
M84 (Rl and R14, respectively). Since each of these re- 
gions is expected to have a contribution from unresolved 
X-ray point sources, a power law component was added 
to the model, which greatly improved the fits (and gave 
a photon index consistent with F w 1.5, as expected, see 
e.g. Sarazin et al. 2003). For regions where this power 
law component could not be well-constrained, we fixed 
the photon index at F = 1.5. Nine regions showed an 
improved fit when another APEC model was added (in 
such cases the abundances of each component were con- 
strained to vary together) . Data from each pointing were 
treated as separate data groups, with normalizations for 
each pointing allowed to vary independently. As an ex- 
ample, the two temperature plus power law fit to the 
spectra for the core of M86 (Rl) is shown in Figure [H] 

Each of the three main components of the diffuse emis- 
sion (the Virgo ICM, M86 group gas, and the cooler cen- 
tral M86 gas) was detectable in some subset of the re- 
gions shown in Figure [71 In R13 it was possible to fit all 
three components simultaneously. The Virgo and M86 
halos cover most, if not all of the FOV. It was there- 
fore necessary to include these components when doing 
detailed spectral fits, even in regions where they could 
not be modeled directly due to inadequate statistics. To 
account for the varying Virgo contribution, we took the 
normalizations of the best fit APEC models to R18 and 
R23 and linearly interpolated across the field based on 
the distance from the center of M87. Although this only 
gives a rough estimate, we note that it does not affect 
our results significantly. For regions closer to M87 (e.g., 
R19) the interpolation gives a value very close to that 
measured from R18, which is an appropriate local back- 
ground. Similarly, regions near the faint end of the M86 
tail (e.g., RIO) have normalizations close to that of R23, 
which again, is an appropriate local background for this 
area. For most regions of interest, M86 is bright enough 
that the contribution from Virgo is negligible, and the 
accuracy of interpolating is therefore sufficient. 

A similar approach was used to model the contribu- 
tion from the halo around M86. Although, in general, the 
normalization of this component decreased with distance 
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from the core of M86, the variation across the field was 
more complicated than for Virgo, which is expected since 
this gas is in the process of being stripped from M86. We 
included a fixed component, with kT = 1.2 keV, abun- 
dance set equal to that of the main component, and vari- 
able normalization, in regions where the M86 halo tem- 
perature could not be determined directly (see column 2 
of Table [2]). The effects of accounting for this component 
can be seen, for example, by comparing lines 1 & 2 in Ta- 
ble [H Including the hotter 1.2 keV component results in 
a lower best-fit temperature for the cooler gas, and gives 
a better statistical fit to the data. In general, our results 
are consistent with Finoguenov et al. (2004), who derive 
X-ray temperatures and abundances for the central re- 
gion, including M86 and the plume from XMM-Newton 
observations. 

We note that for the cores of M86 and M84 (Rl and 
R14 in Table [2]) the abundances are poorly determined, 
despite the large number of net counts. Using a VAPEC 
model to allow for non-solar abundance ratios did not 
improve the fits. The failure of these models may be 
due to abundance variations within the defined regions. 
For instance, the abundance map given by Finoguenov 
et al. (2004) shows significant variation near the core of 
M86, even on relatively small scales. Additionally, the 
models do not allow for multiphase gas, although the 
temperature most likely varies significantly in the cores. 

Using the detailed spectra, we can derive the approx- 
imate mass in cool gas for the base of the tail (the 
"plume") and the core of M86. Using the temperature 
and abundance measured in R4, and assuming a prolate 
spheroid geometry with major and minor axis lengths 
of 12.3 kpc and 7.9 kpc, we find a density of ripiume ~ 
8.2x10"^ cm"^ and a mass Mpiumc ~ 7.4xlO^M0. Sim- 
ilarly, in the core of M86, we find ricoro ~ 6.2 x 10~^ cm"'^ 
and Mcoro ~ 7.4 x 10* Mq within a sphere of radius 
9.6 kpc. Assuming a cylindrical geometry for the tail, 
and considering only the cool gas components in R5 - 
R9, we find a typical density of ritaii ~ 1-5 x 10~^ cm~'^ 
and a total gas mass of Mtaii ~ 1-7 x 10^ M©. The mass 
in stripped gas is therefore more than three times the gas 
mass in the core, which supports the scenario where a sig- 
nificant fraction of the M86 ISM is removed rapidly due 
to ram pressure stripping (Takeda et al. 1984). Rangara- 
jan et al. (1995) find similar results for the plume and 
the core. However, we find consistent abundance values 
in the core and in the plume within the errors, while 
Rangarajan et al. (1995) find a higher abundance in the 
plume. This disagreement is likely due to the multiphase 
model we use, which includes a separate component for 
the 1.2 keV M86 halo emission. 

5. DISCUSSION 

5.1. The Ram Pressure Stripped Tail of M86 

The most striking feature in Figure is the long tail 
of emission extending to the NW from M86. M86 has 
a line-of-sight velocity of -244±5 km s'-^, while M87's 
is 1307±7 km s^^ (Smith et al. 2000). Therefore, M86 
is traversing the Virgo cluster at wm86 > 1550 km s~^ 
(about Mach 2 for kT = 3 keV). The tail naturally forms 
due to ram pressure stripping of the M86 corona by the 
Virgo cluster ICM. Finoguenov et al. (2004) suggested 
that the tail formed due to interactions with a filament 



rather than the Virgo ICM, though their result was based 
on older distance measurements that placed M86 outside 
the virial radius of the Virgo cluster. If we assume that 
M86 is bound to the Virgo cluster, its large line-of-sight 
relative velocity allows us to constrain the separation be- 
tween M86 and M87. Using the M87 mass profile detailed 
in § 15.31 we find that a free-fall velocity of 1500 km s~^ 
corresponds to a separation from M87 of 0.5 Mpc (the 
projected separation is 0.35 Mpc). This is an upper-limit, 
since the total relative velocity must be at least as large 
as the line-of-sight relative velocity. This separation is 
consistent with recent surface brightness fiuctuation dis- 
tances of Virgo cluster galaxies, which give a distance 
between M86 and M87 of 0.4 ±0.8 Mpc (Mei et al. 2007). 

We also can place a lower-limit on the length of the 
long stripped tail using the mass profile of M87 and the 
line-of-sight velocity. The maximum free-fall velocity of 
M86 from infinity at the projected separation of 0.35 Mpc 
is about 1680 km s~^. Using this as an upper-limit on 
its current 3D velocity, we find that the angle between 
the direction of motion of M86 and the line-of-sight is 
9 < 23°. Assuming that the stripped tail is aligned with 
M86's current direction of motion, and given that the 
length of the tail in the plane of the sky is 150 kpc (0.51°), 
we find a lower-limit on the actual length of the tail of 
-^taii ^ 380 kpc, making this the longest ram pressure 
stripped tail presently known. 

The long stripped tail originates in the plume of diffuse 
emission located directly north of M86. For a continu- 
ous stripping process, one would expect the tail to ex- 
tend from M86 itself. This separation between M86 and 
the plume can be explained by rapid recent stripping, 
in which a significant fraction of the remaining gas in 
M86 is rapidly stripped when the ram pressure stripping 
condition is met. A detailed discussion of this process is 
given in § 15.41 

As seen from Figure [7] and Table [21 the stripped tail 
shows a general trend of cooler gas in and near the 
plume (with temperatures in the 0.8-0.85 keV range) 
and warmer gas at the tip (in the 0.9-1.2 keV range). 
The temperature structure of the tail is consistent with 
a ram pressure stripping model, where, as M86 falls into 
the Virgo cluster, the hotter, higher entropy group gas 
is stripped first due to interactions with the Virgo ICM, 
followed by the cooler, lower entropy M86 ISM, which is 
removed rapidly once the stripping condition is met. 

5.2. The Extended M86 Halo 

The Chandra and ROSAT X-ray images (Figure [3]) 
show an extended halo associated with M86, with a sharp 
brightness edge to the southeast in the direction of M87. 
Since M86 is traversing the Virgo cluster supersonically 
(see § 15. ip we expect a shock to be driven in the Virgo 
ICM, possibly producing a brightness edge similar to that 
seen in the southeast. We attempted to measure the den- 
sity jump associated with the shock in the following way. 
We extracted the Chandra 0.5-2.0 keV surface brightness 
profile of this edge in two roughly equal sectors since the 
morphology of this feature is irregular in the south (see 
Figure [3]), and since the gas temperature differs in these 
two regions (see § 14. 1|) . The emission measure profile for 
the northern sector is shown in Figure [9l Distance is 
measured from the center of curvature of the apparent 
edge, which is 10 kpc (2.16') east of M86. This profile 
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was fit with a spherical gas density model consisting of 
two power laws. The free parameters were the normal- 
ization, the inner (a) and outer (/?) slopes, the position 
of the density discontinuity (rbroak), and the amplitude 
of the jump {A). We assumed that the gas is isothermal 
with kT — 1 keV and that the abundance is constant 
at 30% solar, consistent with results from spectral fits 
(see § [4]) . For the best fit model (see Figure ^ , we find 
a = 0.491^:1^, 13 = -0.8lt°;i^, rbrcak = 521^ kpc, and 
A — 1.3to;4- Similar results are found from the southern 
sector, but with larger errors. The lack of a well-defined 
edge is consistent with our findings for the orbit of M86 
(see ij 15. 3p . In particular, for all of the likely orbits, M86 
is moving to the southeast, such that our lines of sight 
pass through the Mach cone. To see a sharp edge, our 
line of sight must be tangent to the shock front. 

5.3. Constraints on the Orbit of M86 

Knowing the orbit of M86 is the key to understanding 
its interaction with the Virgo cluster. The ram pres- 
sure stripped tail reveals the motion of the galaxy on the 
sky. The length and direction of the tail, together with 
the large line-of-sight speed of M86, constrain its orbit. 
The large luminosity of M86 suggests that it dominates 
the associated in-falhng subgroup (Schindler et al. 1999). 
In the following, we assume M86 is bound to the Virgo 
cluster. For the purpose of calculation, the gravitational 
potential of the Virgo cluster is treated as a spherical 
NFW potential (Navarro, Frenk & White 1997), with 
a virial radius of 1.3 Mpc (Evrard, Metzler & Navarro 
1996) and a concentration parameter of 4.5 (Neto et al. 
2007), appropriate for a ~ 3 keV cluster. The gravitat- 
ing mass was normalized to match the total mass within 
320 kpc of 4.4 X lO^^ (Schindler et al. 1999; scaled 
to a distance of 16 Mpc). For this potential, the escape 
speed from a radius of 351 kpc (the projected separa- 
tion of M86 from M87) is 1677 km s~^, not much larger 
than M86's line-of-sight speed of 1550 km s^^ relative 
to M87. The empirical mass distribution of Schindler 
(1999) gives the same escape speed at 351 kpc as our 
model potential. We note that the gravitational poten- 
tial of the Virgo cluster is not well constrained at dis- 
tances from the cluster center comparable to or greater 
than its virial radius. Furthermore, the truncated NFW 
potential is highly simplified, ignoring mass beyond the 
virial radius, departures from spherical symmetry, and 
the dynamic state of the cluster due to continuing in- 
fall. Since some of the limits we derive here are sensi- 
tive to the poorly known potential at large radii, they 
should be treated as indicative rather than quantitative 
(outside the context of the model). The material in this 
section is supplemented by a more detailed discussion in 
Appendix [XI 

We consider radial orbits first, as suggested by previ- 
ous studies (e.g., Forman et al. 1979; White et al. 1991). 
In order for the gas tail to point away from the cluster 
center, the galaxy must be inbound. It can easily be 
shown that the line-of sight velocity, fios, is maximized 
for some radius r greater than the observed separation 
s (see Appendix \Q . For a marginally bound (zero en- 
ergy) radial orbit in the potential described above, the 
maximum value of |wios| occurs when M86 is close to 
twice its projected distance from the cluster center, giv- 



ing |t;ios| — 1205 km s^^, less than the observed value of 
1550 km s~^. Thus, unless M86 is significantly unbound 
from the Virgo cluster (or the potential is incorrect), its 
line-of-sight speed is inconsistent with radial and nearly 
radial orbits. 

It is convenient to specify more general orbits in terms 
of their inner and outer turning radii, rj and Tq, respec- 
tively. We constrain the possible orbits of M86 by placing 
limits on these parameters. Coarse limits can be placed 
on these radii by considering energy arguments alone. 
Assuming that M86 is bound to the Virgo cluster, we 
find n < 489 kpc and Tq > 3.5 Mpc, or ~ 2.7 times 
the virial radius of the Virgo cluster (see Appendix [A|). 
These limits can be further restricted by considering the 
range of possible viewing directions for each point on an 
orbit. There can be zero, two, or four possible viewing 
directions that would place M86 at the observed separa- 
tion from the cluster and give it the observed line-of-sight 
velocity (see Appendix \K\ . Figure [TO] shows a range of 
marginally bound orbits [r^ ~ oo), with regions that 
meet these conditions marked in color. The range of an 
orbit where these conditions are met shrinks as Tq de- 
creases, i.e., as the orbit becomes more tightly bound. 
For the marginally bound orbits, the full range of for 
which v\os can attain its observed value is 247 < rj < 395 
kpc (the range represented in Figure llOp . For more 
tightly bound orbits, the acceptable range of rj is re- 
duced. 

The lower limit on is increased if we consider the ex- 
tent of the ram pressure stripped tail. The projected or- 
bit must extend to at least the distance (~ 100 kpc) that 
the tail projects beyond M86 in the direction away from 
the cluster center. Locations on the orbits where the line- 
of-sight speed can attain —1550 km s^^ and these con- 
ditions also are met are shown in red in Figure 1101 We 
see that the range of possible locations for M86 on these 
orbits is tightly constrained. We repeated this analysis 
for a number of values for the outer turning radius. The 
range of potential orbits and locations for M86 shrinks 
with decreasing and no orbits were found to meet these 
criteria for Tq < 8.2 Mpc (6.3 virial radii). 

Figure [Tl] shows projections onto the sky of the 
marginally bound orbits (corresponding to the midpoints 
of the red regions of Figure [TO]) overlaid on the Chandra 
0.5 - 2 keV image, with M86 at the representative loca- 
tion for each orbit. It is evident that the possible orbits 
for M86 can be pruned further. For example, the or- 
bit on the lower right in Figure [Tl] does not reach far 
enough north to produce the remote part of the gas tail. 
At the other limit, the position angle of M87 measured 
from M86 is 116° (east from north), generally eastward 
of these orbits. Overdense stripped gas tends to fall to- 
wards the cluster center, in that direction, so that some 
stripped gas can reasonably lie to the east of the orbit. 
Stripped gas cannot lie west of the orbit. Thus the orbits 
that cross the region to the east of the tail and north of 
M86, where there is no sign of stripped gas, are unlikely 
candidates for the orbit of M86. 

Note that the orbits in Figure [TT] are not simply related 
to their inner turning radii. The orbit corresponding to 
the smallest value of the inner turning radius appears 
second from the right at the top of the figure. At first, 
orbits for increasing values of n lie to the left of this, but 
at a value of r^ approaching (but less than) s, the accept- 



6 



RANDALL ET AL. 



able viewing direction flips from outward to inward (i.e., 
the acceptable sign of a' changes, see Appendix [A|) . At 
this point, the projected position of the orbit shifts from 
leftmost to rightmost at the top of Figure [TlJ Whereas 
the apparent location where the orbit passes through the 
virial radius was moving outward, it now moves inward as 
Ti increases, coming into the field of view for the largest 
values of ri here. This results in disjoint ranges of pos- 
sible orbits for M86. Thus the marginally bound orbits 
with inner turning radii close to rj ~ 260 kpc and those 
with inner turning radii in the range 330 "^r-i < 380 kpc 
are best suited to model the orbit of M86. 

Applying the same criteria to the more tightly bound 
orbits yields the collection of possible orbits shown in 
Figure fT2l We have not attempted to be highly selective 
or exhaustive, since that requires a well defined model 
for the formation of the gas trail. For values of the outer 
turning radius in the lower end of the acceptable range, 
the shrinking range of possible orbits and locations ex- 
cludes orbits with inner turning radii rj > s. All of the 
acceptable orbits for M86 are weakly bound to the Virgo 
cluster, with the most tightly bound having outer turn- 
ing radii at 8.8 Mpc. All of the possible locations for 
M86 he only a little farther from M87 than M86 does in 
projection. All are close to the plane of the sky, ranging 
from 167 kpc closer than M87 to 263 kpc farther than 
M87 from the Sun at the extremes. M86 must also be 
close to the pericenter of its orbit. The direction of mo- 
tion of M86 is close to our line of sight, within 16° - 23° 
of it, for all of the cases illustrated. 

M86 is traversing the Virgo cluster supersonically 
(I'los = —1550 km s~^ alone is almost twice the sound 
speed Cs, which is taken to be 850 km s~^). Therefore, we 
expect a shock to be driven in the Virgo ICM. For the or- 
bits we derive, the angle of inclination is smaller than the 
opening angle of the Mach cone (~ 33°) and we should 
not expect shock fronts to be visible on the plane of the 
sky (cf. Rangarajan et al. 1995). This is illustrated in 
Figure fT3l which shows the marginally bound orbit with 
Ti = 376 kpc. Each circle is drawn centered on a point 
where M86 was at a time 5t in the past, with a radius of 
Cs 6t. This is done at equally spaced times to indicate the 
shape of the Mach cone. The absence of caustics in Fig- 
ure[T3]shows that we would not see the shock front in pro- 
jection. (While this statement is accurate, except possi- 
bly in the small region where the shock front is highly su- 
personic, the diagram is schematic. For example, if M86 
were coming directly towards the Earth, radii of sections 
of the Mach cone would be (1 — Cg/u^)^^/^ ~ 1.20 larger 
than drawn.) Our lines of sight through the compressed 
gas behind the shock are longest where the circles pile 
up to the southeast of M86, consistent with the enhance- 
ment in X-ray surface brightness seen ~ 3' southeast of 
M86 (Figure [131). 

5.4. Displacement of Stripped Gas 

The prominent plume of gas lying 3'- 4'(~ 16 kpc) 
north of M86, noted previously (e.g., Forman et al. 
1979), appears to be physically separated from the dense 
gas remaining in the galaxy and it is centered some dis- 
tance east of the orbits illustrated in Figure [T^ The re- 
moval of a significant fraction of the ISM in a single blob 
is expected from rapid ram pressure stripping (Takeda 
et al. 1984). A highly simplified model for ram pressure 



stripping treats the plume as a single particle subject 
to gravity and drag due to its motion through the sur- 
rounding gas. Allowing for buoyancy, the net force on the 
plume due to gravity is {p — pc)V g„ where p is its density, 
V is its volume, pc is the density of the ambient ICM, 
and g is the local acceleration due to gravity. The drag 
force on the plume is — Cd^PoWV, where Cd is the drag 
coefficient, A is the plume cross section, v is its velocity, 
and V = |v|. Thus, the net acceleration of the plume is 
a = [1 — Pc//o]g — [CtiPcA/{pV)]vv. Two of the main 
parameters of this model are the density contrast, pd p, 
and the factor C-qA/V . For the density of the ICM, po, 
we use the beta model of Schindler et al. (1999; /3 = 0.47, 
core radius ~ 2.7'). On the grounds that the density of 
the plume is likely to be lower now (due to its stripping 
and ejection from the confining potential of M86), we set 
its density contrast at the projected radius of M86 to be 
18 (cf. ~ 16 from the numbers above, see ij |4.2|) and treat 
the density of the plume as constant. The other factor 
is determined by treating the plume as a sphere with a 
constant radius of 10 kpc, with Cd = 0.75. As above, 
this radius is a little larger than the observed radius. To 
simulate stripping, a gas blob is placed at the location 
and velocity of M86 at the time it is released from the 
cluster virial radius. The orbit of the blob is followed, 
subject to the gravitational acceleration of the cluster 
and the moving galaxy. 

This model cannot account for the current location of 
the plume if the gravitational potential of M86 is spher- 
ical. The path of the blob is determined by the compe- 
tition between gravity and drag. If the drag is large, the 
blob is ejected from M86 early on its orbit, and slows 
quickly until it is falling towards the cluster center at 
its terminal speed, and ultimately lies farther away from 
M86 than observed. Its early ejection also leaves it too 
far back along the orbit of M86. Reducing the relative 
significance of the drag causes the blob to be ejected 
later, bringing it closer to its observed position but not 
far enough from the orbit in the direction of the cluster 
center. 

This issue can be resolved if the gravitational potential 
of M86 is aspherical. Consider a small gas blob in an as- 
pherical galaxy moving at an inclined angle through the 
ICM, illustrated schematically in Figure [HI While ram 
pressure is insufficient to eject the blob, it is driven to 
an equilibrium position where ram pressure is balanced 
by gravity. Since the direction of the gravitational field 
must oppose the drag, the external flow must be per- 
pendicular to the equipotential surface. The equilibrium 
position therefore lies at a location away from the axis 
passing through the center of the galaxy and parallel to 
the flow (e.g., near the point labeled "equipotential" in 
Figure [T4[) . As the drag starts to overwhelm gravity and 
displaces the blob in the direction of the flow, the equipo- 
tentials the blob encounters tilt, so that it is subject to a 
component of the gravitational force directed away from 
the axis. Because forces in the direction of the flow are 
nearly balanced, this off-axis component of gravity drives 
the blob farther away from the axis of the flow. Thus, a 
blob stripped from such a galaxy tends to emerge away 
from the axis of the flow. Treating the gas as a fluid, 
interstellar gas that is pushed inward by ram pressure 
along the leading edge of the galaxy flows out preferen- 
tially along the long axis of the potential, in the direction 
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of the weakest gravitational force. 

The model for the aspherical gravitational potential 
of M86 used here has the form $(r) = F{'w), where 
F(w) = '01n[(l + w/c)/{\ + w)\/w is a minor modifi- 
cation of the standard NFW form (■!/'ln[l/(l + w)]/'w in 
the same notation). The extra factor of 1 + w/c, where 
c is the concentration parameter, makes the total mass 
converge, avoiding the need to truncate the mass dis- 
tribution at the virial radius (which would create prob- 
lems for our fourth order integrator). The coordinate 
w = \J [x^ + y'^)/a? -\- l\P- + 1, where a is the NFW 
potential scale length and h is chosen to give the desired 
ellipticity. The additional 1 under the square root flat- 
tens the potential at small r, providing a better model 
for the gravitational force on an extended gas blob (of 
size comparable to a) when it is close to the center of 
the galaxy. This modification also makes the integrator 
behave better near the center of the galaxy. The virial 
radius of M86 was set to 100 kpc and its concentration 
parameter to 8, roughly the values expected for a mas- 
sive galaxy. The normalizing factor, -0, was expressed as 

= 9(7^, so that CT is a rough measure of the line-of- 
sight velocity dispersion for M86. The value used below 
is CT — 185 km s~^. We note that for the model what 
matters is the ratio of the drag force on the blob to the 
binding force of the M86 potential, allowing for a trade 
off between the choice of potential and the density ratio 
and size of the blob. As a result, the choice of potential 
is not very critical. 

The two remaining parameters are the axial ratio, a/5, 
and the orientation of M86. Consistent with its S0/E3 
classification, optical images show that M86 is highly flat- 
tened on scales comparable to the projected distance of 
the plume (e.g. a/6 ~ 1.56 for the isophote with a ~ 4.5' 
in the deep image of Nulsen & Carter 1987). Based on 
the discussion above, the potential of M86 would need 
to have a long axis pointing roughly northward from our 
line-of-sight. Therefore, we assume it is oblate, using 
a/5 = 2. The orientation of M86 is then determined by 
the direction of its minor axis. From the optical images, 
this has a position angle of ~ 35° on the sky. However, 
its tilt with respect to the plane of the sky is unknown. 
We have set the minor axis to point towards us in the 
north, at an angle of 45° from the plane of the sky. This 
is roughly the orientation that maximizes the transverse 
displacement of the blob. Figure [15] shows the path of the 
blob on the sky for this set of model parameters for the 
orbit with Tq = 9.1 Mpc and n = 314 kpc. In this model, 
the blob is currently slower than M86 by 355 km s~^ 
along our line of sight and trails it by 38 kpc. 

Around the edges of the gas halo in M86, where the 
shear in the external flow is strong (Figure [M]) , if the 
effective viscosity is high, viscous stresses pull the inter- 
stellar gas out of the galaxy. Alternatively, if the vis- 
cosity is low, shear instabilities mix the interstellar gas 
with the ICM, also stripping it from the galaxy (Nulsen 
1982). This stripping is aided by the low pressure due to 
the Bernoulli effect around the edges of the galaxy, which 
tends to pull gas into the path of the flow. This process 
works around all edges of the inclined galaxy, as seen 
from the direction of the flow. However, it is expected to 
be greatest at the leading and trailing edges. The shear 
in the external flow is expected to be greatest near the 



leading edge, favoring stripping there. The large ram 
pressure at the leading edge displaces the edge of the 
interstellar gas deeper into the potential of the galaxy, 
pushing it towards the trailing edge of the galaxy, where 
the pressure is lower. Thus, gas at the trailing edge sits 
higher in the gravitational potential of the galaxy, favor- 
ing its removal from the galaxy. This may account for 
the apparent double streams of gas seen trailing M86 and 
several other galaxies in the composite image. "Viscous" 
stripping from the main body of M86 as well as the plume 
can explain the broad features of the gas trail, though 
modeling the finer features, such as the smaller blobs 
of gas lying to the east of the galaxy, requires a more 
detailed treatment of the gas dynamics (e.g., numerical 
simulations) and is beyond the scope of this paper. 

The proton mean free path due to Coulomb collisions 
is given approximately by A ~ 150(fcT)^n~3 pc, where 
kT is the gas temperature in keV and the electron den- 
sity is 10~^n_3 cm"'^. In terms of this, the Reynolds 
number is Re ~ wL/(sA), where v is the flow speed, s is 
the sound speed, and L is a relevant length scale. Tak- 
ing V — 1550 km s^^ and L — kpc gives Re ~ 20 
in the 2.4 keV ICM, but Re ~ 3000 for the gas in the 
0.77 keV plume (and interstellar gas) if it is exposed di- 
rectly to the external flow. These values suggest that 
the external flow can be largely laminar, while the flow 
in the cool interstellar gas is relatively turbulent. Fur- 
ther complicating matters, the effective mean free path 
may be significantly smaller than the Coulomb mean free 
path (e.g., Schckochinin et al. 2007). The Reynolds num- 
ber for this flow is therefore not well defined, due to the 
range in gas temperature and uncertainty in the mean 
free path. 

6. SUMMARY 

We have argued that the plume and long tail of M86 
formed due to ram pressure stripping forces generated 
as M86 falls into the Virgo cluster. Several studies have 
found a similar interpretation for the formation of this 
feature (e.g., Forman et al. 1979; Fabian et al. 1980; Ran- 
garajan et al. 1995, however see Bregman & Roberts 
1990; Finoguenov et al. 2004). We concentrate on these 
main results: 

• the plume and long tail observed in the diffuse 
emission are created by ram pressure stripping as 
M86 falls into the Virgo cluster. The tail is 150 kpc 
in projection (a simple estimate, which assumes 
free-fall velocity for M86 and an NFW potential 
for M87, gives a lower-limt on the true length of 
the tail of 380 kpc), making this the longest ram 
pressure stripped tail presently known. 

• based on the X-ray spectra, we detect three distinct 
components associated with the M86/Virgo clus- 
ter system: the Virgo ICM, with kT ^ 2.4 keV; 
the extended halo of M86, with kT ^ 1.2 keV; 
and the cooler central and stripped gas of M86, 
with kT ~ 0.8 keV. The temperature structure of 
the tail is consistent with ram pressure stripping, 
where the higher entropy M86 halo gas is stripped 
first and deposited in the tip of the tail, and the 
lower entropy M86 ISM is stripped more recently, 
constituting the base of the tail and the plume. 
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• the large line-of-sight velocity of M86, and its posi- 
tion relative to the Virgo cluster, tightly constrain 
its orbit, especially if it is assumed that the gas 
tail traces the orbit. In particular, the observa- 
tions are inconsistent with a radial orbit. We show 
that M86 is at best only marginally bound to the 
Virgo cluster, with an inner turning radius on the 
order of 300 kpc as expected from its recent in-fall. 
Our best-fitting orbital model requires that M86 
be close to M87, less than 167 kpc closer than or 
263 kpc farther than M87 along our line of sight, 
which is consistent with the most recent distance 
estimates based on surface brightness fluctuations 
(Mei et al. 2007) which give a line-of-sight separa- 
tion of 0.4 ± 0.8 Mpc. 

• the prominent plume of gas lying 3'- 4'north of 
M86 appears to have been rapidly driven from M86 
by ram pressure stripping. The projected position 
of the plume, which does not lie directly on our 
best-fit model orbit for M86, can be understood 
if M86 has an aspherical potential (as suggested 
by optical isophotes). If M86 moves through the 
Virgo ICM at an inclination angle relative to the 
local flow, the gas at the trailing edge is more eas- 



ily stripped, thereby displacing the gas from the 
nominal orbit of M86 itself. This model may also 
explain the apparent double streams of gas seen 
trailing M86, as well as those in other Virgo galax- 
ies. 

• the apparent brightness edge to the southeast seen 
in i? 05"^ T" observations is also seen in the Chandra 
images. The edge is well fit with a two power law 
gas density model, with an abrupt jump in density 
by a factor of 1.3lg;4 edge (consistent with 

no jump). Assuming that this brightness edge is 
the shock generated as M86 supersonically falls into 
the Virgo cluster, the lack of a well-defined density 
jump is consistent with what is expected from our 
orbital model, which suggests that the orientation 
of the Mach cone would make it difficult to detect 
the shock edge. 
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TABLE 1 
Observation Details 



yjos lu 


uate vJDs 


larget 


rt.cTive o^us 


a 

Oleaned Exposure 
(ksec) 


318 


2000-04-07 


M86 


S3, S2, 12, 13 


12.964 


803 


2000-05-19 


M84 


S3, S2, 12, 13 


26.699 


963 


2000-04-07 


M86 


S3, SI, S2, 12, 13 


13.220 


1619 


2001-06-08 


NGC4388 


S3, S2, 12, 13 


19.705 


2882 


2002-01-29 


NGC4438 


S3, S2, 12, 13 


24.891 


5908 


2005-05-01 


M84 


S3, SI, S2, 12, 13 


36.299 


5912 


2005-03-09 


SE of M86 


10, 11, 12, 13, S2 


31.370 


5913 


2005-03-19 


M86 Tail Tip 


10, 11, 12, 13, S2 


34.744 


6131 


2005-11-07 


M84 


S3, SI, S2, 12, 13 


37.855 



* Total cleaned exposure is ~238 ksec. 
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TABLE 2 
Spectral Fits" 



Region # 


kT 


Abund. 


r 


9/1 c 

X /dof 


Net Cnts. 


{Key ) 








1 


0.724i;!;0i2//(2.4) 




-'-■'-'--0.68 


108/91=1.18 


6235 


1 


0.677lH^5/1.29lti;;«S/(2.4) 




(1.5) 


90/89=1.01 


6235 


2 


0-697t»;;!;5S/(1.2)/(2.4) 






56/57=0.98 


3086 


3 


0-653±E:::;«/(1.2)/(2.4) 




(1.5) 


88/83=1.06 


4273 


4 


0.773lf,;^i|/{1.2)/(2.4) 


-0.12 




109/84=1.29 


8630 


5 


0-854t«;'«6/(l.2)/{2.4) 


'-'■^'J-0.09 




94/93=1.01 


4461 


6 


0.849t«;;!S/(i.2)/(2.4) 


'J-'^J^-o.io 




89/88=1.01 


5059 


7 


o-86o±;;;;!«j/(i.2)/(2.4) 


'^■^■^-0.20 




87/74=1.18 


2843 


8 


0-785l»;5;^^/(1.2)/(2.4) 


"^•^3-0. 17 




44/40=1.09 


1749 


9 


0-924t°i|i/{1.2)/(2.4) 


o.solf,:™ 




30/29=1.03 


1043 


10 


1.19liO;J27/(2.4)b 


^^-0.25 




20/25=0.73 


739 


11 


o.865i;;:i»'^/i.29oii;:?«/(2.4) 


'J-'^'^-O.l,'-) 




163/176=0.92 


8795 


12 


0.87ll«:J»;?/1.176l»:?g/(2.4) 


n 62+"''^*^ 

'^■"^-0.19 




240/187=1.29 


8914 


13 


675+°'"l /I 256+°'^^'' /2 723+0-7*-* 


n 72+"-48 

"•'^-0.18 


1 09+0.11 


343/293=1.17 


11033 


14 


0.625l«:««/(1.2)/(2.4) 
0.790i;;-'i';;i!/(1.2)/(2.4) 


1 17+0.78 
-0.10 


459/327=1.40 


42002 


15 






93/90=1.03 


3666 


16 


1.355t»:;^«/(2.4)'> 


0.49+E!-?? 




57/59=0.96 


1746 


17 


1.09lt»:?S/2.422lf,:g«'> 


'J-^'J-0.20 




53/45=1.18 


3054 


18 
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Each region contained up to three thermal APEC components: a low temperature, a mid-range temperature, and a high temperature. 
Temperatures are given, in that order, in column 2. The high temperature component models background Virgo ICM emission. Where 
this component could not be accurately measured its temperature was fixed at 2.4 kcV (sec text for discussion). Similarly, the mid-range 
temperature component models M86 halo emission, and where the temperature could not be measured it was fixed at 1.2 kcV. The 
abundance for this component was tied to the abundance of the cooler fitted component. A power law component was included when 
necessary, presumably for unresolved point sources. Galactic absorption was assumed throughout.'' No low temperature component was 
included."^ Upper bound of 90% confidence interval could not be determined. 
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Fig. 1. — Exposure corrected, background subtracted 0.5-2 keV mosaic image of the Chandra observations in the region of M86. The 
image has been smoothed with a 6" radius gaussian. For each pointing, regions with less than 10% of the total exposure for that observation 
were omitted. Complex structure in the diffuse emission is observed near and in the stripped tail of M86. Directly to the west of M86, 
M84 is visible, which, aside form showing complicated structure in the core, shows a tail of diffuse emission to the south. Also visible are 
NGC 4388, to the south of M86 and M84, and NGC 4438, directly east of M86, both of which show structure in their diffuse emission. 
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Fig. 2. — Same as Figure [T] but with different scaling and a larger smoothing kernel (9") to highlight the structure at the low surface 
brightness end of the M86 tail and to show the faint M84 tail more clearly. Point sources have been removed (see text for details). The 
main M86 tail, after turning from northwest to north, appears to turn back to the west near the very tip. 




Fig. 3. — (left) ROSAT PSPC view of M86. Diffuse emission from M87/Virgo is clearly visible in the southeast, as is an extended halo 
of diffuse emission associated with M86. Also shown are the Chandra (center) and DSS (right) images of the same region. 
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Fig. 4. — Close up view of the long stripped tail, in the 0.5-2 kcV band and smoothed with a 16" radius gaussian. The tail clearly turns 
from a NW direction to the north, and back directly east at the very tip. Furthermore, there appears to be a faint secondary tail along 
the northern edge of the faint end of the main tail, which traces its path. The red box indicates the region used to generate the plot in 
Figure \5\ 
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Fig. 5. — Plot of the total count rate in evenly-spaced binned regions taken across the width of the box shown in Figure |4] The x-axis 
gives the R.A. of each bin on the southern edge of the box. 
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Fig. 6. — Left Panel: A closer view of the central region of Figure [2] with scaling chosen to highlight structure in the central region 
and nearby diffuse emission. The region used to extract the temperature map is overlain. Center Panel: Derived temperature map of the 
same region. The contribution from the Virgo ICM has been taken out. Zero- value d pixels indicate areas that required regions larger than 
the selected maximum area in order to include the minimum net counts (see ^ 14. 11 1. Pixels in faint regions which gave an anomalous fitted 
temperature were also removed. The colorbar under the right panel gives the temperature in keV. Right Panel: Tessellated temperature 
map. Each bin contains roughly 1100 not counts. The counts in each bin were fit with a single APEC model. The contribution from the 
Virgo ICM has not been taken out, leading to higher temperatures farther from M86 as compared to the center panel. The size of the bins 
roughly indicate the size of the extraction region for each pixel in the center panel. 
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Fig. 7. — Same as Figure |2l but with the regions used for extracting spectra overlain. Single regions for spectral fitting are shown in 
green, while the regions shown in white arc fit as one combined region. Table [2] gives the region number and the best-fit tempcrature{s). 
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Fig. 8. — Data (crosses) and fitted model (histogram) for the core of M86, Rl in Table [2] The model is a two-temperature thermal 
plasma with fitted temperatures and abundance (given in the second row of Table [2j, plus an unresolved power law component with photon 
index of 1.5 and a fixed thermal plasma component to model emission from the Virgo ICM, all with Galactic absorption. The data are 
from two S3 observations, OBS-IDs 318 and 963 (see Table[TJ. 
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Fig. 9. — Emission measure profile of the northern half of the bright edge seen in the extended group omission southeast of M86 (shown 
most clearly in the ROSAT image in the left panel of Figure [Sjl- The x-axis gives the radius from the apparent center of curvature defined 
by the feature. The best fit two power law density jump model is given by the solid line. 
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Fig. 10. — Possible marginally bound orbits for M86. Orbits are plotted for inner turning radii in the range 247 - 395 kpc. The direction 
of motion is right to left. The dotted circle shows points 351 kpc from the cluster center. The colored regions (dash-dot and solid lines) 
show where on each orbit the galaxy could be viewed at a projected distance of 351 kpc from the cluster center, with a line-of-sight speed 
of 1550 km s~^ towards us. Locations on the orbits where the past orbit would also appear to project 100 kpc outside the current position 
of M86 are shown in red (solid lines; the midpoint of each of these regions is marked with an asterisk). Both regions shrink for orbits with 
smaller outer turning radii. As plotted, all orbits pass through the virial radius at x = 1300 kpc on the +x axis. 
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Fig. 12. — Possible orbits for M86. As in Figure [TTI one representative orbit is shown projected onto the smoothed 0.5 - 2 keV image of 
M86 for the remaining plausible orbits. Colors indicate the outer turning radius as follows: To = 8.76 Mpc, yellow; To = 9.11 Mpc, blue; 
To = 19.3 Mpc, green; To = oo, red. 
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Fig. 13. — Illustration of the Mach cone for M86. Shown in yellow is the marginally bound orbit with ri = 376 kpc, one of the better 
candidates for the orbit of M86. Centered on the orbit, at equally spaced times, circles are plotted on the plane of the sky with radii Cs St, 
where St is the time before the present when M86 was at that location and Cg is the sound speed (850 km s~^). 
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Fig. 14. — Stripping from an aspherical galaxy. Schematic diagram of an inclined, aspherical galaxy moving supersonically to the left 
through an external medium. The solid ellipse represents an equipotential of the galaxy, the dotted line represents the shock front and 
the dashed line represents the outer boundary of the interstellar gas. The ram pressure is maximized near the leading edge of the galaxy, 
where the shock is normal to the external gas flow, and lower on inclined sections of the shock front. The gas flow is fastest around the 
edges of the galaxy, causing local minima in pressure that pull gas out sideways from the galaxy. Gas displaced by ram pressure from the 
leading edge of the galaxy piles up along the trailing side and edge. 
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Fig. 15. — Model orbit for ram pressure stripped blob. The orbit of M86 is shown in white and the (coarsely sampled) orbit of the blob 
is shown in yellow. Note that the path of the blob follows a trail of denser gas. 
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APPENDIX 

CONSTRAINING M86'S ORBIT 

In cylindrical polar coordinates, with the origin at the cluster center and the orbit (instantaneously) in the x - y 
plane, the position of M86 can be expressed as r = re^ and its velocity as v = u^e^ + v^e^, where e^, e^, are the 
coordinate basis vectors. When M86 is viewed along the unit vector, n, its location projected onto the plane of the 
sky is r — n • rn. Expressing n in terms of the coordinate vectors as n = o'e^ + b'e^ + c'ez, the projected separation, 

s, of M86 from the cluster center is given by = (r — n • rn)^ = r^(l — a'^), which requires a' — ±-y/l — s^/r^ (and 
r > s). Since n is a unit vector, a'^ + b'^ + c'^ — 1, further requiring b' ~ i^s^/r^ — c'^. The remaining parameter, 
c', is the cosine of the polar angle of the viewing direction (i.e., the orbit inclination). With the two values of a' above, 
the parametrization {b',c') = (cos ^, sin ^)s/r, for < ^ < 2tt, specifies the complete range of viewing directions that 
project MSG at the observed distance from the cluster center. In this notation, the velocity of M86 along our linc-of 
sight is wios = n • V = a'vr + b'v^. For radial orbits with M86 inbound (vr < 0), wios = = -^/l — jr'^Vr < 0. 
This vanishes for r —>■ s (where a radial orbit must be perpendicular to our line-of-sight), while the radial speed, \vr\, 
decreases with increasing r, so that |wios| is maximized for some r > s. 

The inner and outer turning radii, n and Tq, respectively, are related to the specific energy of an orbit hy E = 
(r^^o — rf(j)i)/{rf — r^), where (pi = (/>(n) is the gravitational potential at n and (po is that at Tq. The specific angular 
momentum of the orbit, £, is given by i'^ = 2rfr1{4)o — 4>i)/{^o — i^f)- If M86 is bound to the Virgo cluster, then, 
since r > n and its speed, v > |wios|, > E — w^/2 + 0(r) > + (pi, which requires (f>i < Taking 

uios — —1550 km s^^, this gives an upper limit on the inner turning radius of n < 489 kpc. To obtain a lower limit on 
ro, note that dE/dvo — (w^ ^ — v1)ro/ (r^ — r?), where the Kepler speed at is given hy v"^^— r d(f>/dr\^^^ and the 

speed at the outer turning radius is = i'^/r'^ = 2{E — (j)o)- The net acceleration at the outer turning radius is inward, 
requiring Vo < wk,o, so that E is an increasing function of Tq. Interchanging and in this argument and noting 
that the net acceleration is outward at the inner turning radius (wK,i < ^^i) shows that E is also an increasing function 
of rj. Since dE = dE/dridri + dE/dradvo, for a fixed energy, maximizing n minimizes Tq. For M86, v > \vios\ and 
r > s, so that its energy E > E^nin = ^los/^ + 4>{s)- To obtain E = E^min with ri = s requires — ro.min — 3.5 Mpc, 
or ~ 2.7 times the virial radius of the cluster. For values of rj larger than s, the total energy would need to match 
or exceed vf^^/2 + (pi > i?min- The value of required to attain this energy is an increasing function of ri, unless the 
mean density of the cluster increases with radius. Thus, ro^min provides a lower limit on the outer turning radius of 
M86's orbit. 

At a fixed point on an orbit, the requirement vios = a'vr + b'v^ determines the possible viewing directions, if any, 
that could model M86. For each of the two solutions a' — zb^/l — s^/r'^, the requirement Vios = a'vr + b'v^ gives 
b' — {vios — a'vr)/v^. If either of these gives b' in the range [—s/r,s/r], then it corresponds to two possible viewing 
directions, with c' — zt-y/s^/r^ — 6'^ (symmetrically placed above and below the plane of the orbit). Thus, for each point 
on an orbit, there can be zero, two, or four potential viewing directions that would place M86 at the observed separation 
from the cluster center and give it the observed line-of-sight velocity. Figure[TO]shows a range of marginally bound orbits 
(ro = cxd), with regions that meet these conditions marked in color. Since |fios| < wios.max — \J'^ — s^/r^jwr -I- s/rjw^l, 
which is an increasing function of for fixed r and rj, the range of an orbit where these conditions are met shrinks 
as To decreases, i.e., as the orbit becomes more tightly bound. 

If the trail of stripped gas lies along the orbit of M86, then the past orbit must project away from the cluster 
center at least as far as the gas tail extends. From above, (r — n • rn)/s is a unit vector in the plane of the sky 
that points from the cluster center towards M86. For each potential location of M86 on an orbit, we have computed 
(i(r') = r' • (r — n • rn)/s for values of r' ranging over the past orbit. If stripping is to account for the location of 
the gas tail, (i(r') must increase initially as r' moves backward along the orbit and the maximum value of (i(r') — s 
must exceed the distance, ~ 100 kpc, that the tail projects beyond M86 in the direction away from the cluster center. 
Locations on the orbits where the line-of-sight speed can attain —1550 km s~^ and these conditions are also met are 
shown in red in Figure [TUl 



